Gas-kinetic-based traffic model explaining observed hysteretic phase transition 
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Recently, hysteretic transitions to 'synchronized traffic' with high values of both density and traffic 
flow were observed on German freeways [B. S. Kerner and H. Rehborn, Phys. Rev. Lett. 79, 4030 
(1997)]. We propose a macroscopic traffic model based on a gas-kinetic approach that can explain 
this phase transition. The results suggest a general mechanism for the formation of probably the 
most common form of congested traffic. 
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To physicists, non-equilibrium phase transitions are 
very fascinating phenomena. Prominent examples are 
pattern-forming transitions in hydrodynamic systems 
driven far from equilibrium, like thermal convection of a 
fluid heated from below or transitions to a state of spatio- 
temporal chaos [0 . Recently, physicists got interested in 
the spatio-temporal, collective patterns of motion formed 
in social or biological systems of so-called 'motorized' or 
'self-driven' particles gj. A particularly strong physi- 
cal activity has developed in the rapidly growing field 
of traffic dynamics P-|l6|], not only because of the large 
potential for industrial applications. 

On a macroscopic scale, many aspects of traffic flow 
are similar to those of aggregated physical systems. In 
particular, if one abstracts from the motion of the single 
vehicles, traffic can be modelled as a continuum com- 
pressible fluid QJ^] (see Ref. M for an overview). Exist- 
ing macroscopic traffic models have been able to explain 
various empirically observed properties of traffic dynam- 
ics, including the transition of slightly disturbed traffic 
to traffic jams ('local cluster effect') [[7j. 

Recently, Kerner and Rehborn presented experimental 
data indicating a first-order transition to 'synchronized' 
traffic (ST) Traffic data from several freeways in Ger- 
many HU and the Netherlands indicate that ST is 
the most common form of congested traffic. ST typically 
occurs at on-ramps when vehicles are added to already 
busy 'freeways' and has the following properties: (i) The 
dynamics of the average velocities on all lanes is highly 
correlated ('synchronized'), (ii) ST is characterized by a 
low average velocity, but, in contrast to traffic jams, the 
associated traffic flow is rather high, (iii) The transition 
to ST is usually caused by a localized and short pertur- 
bation of traffic flow that starts downstream of the on- 
ramp and propagates upstream with a velocity of about 
— lOkm/h. (iv) As soon as the perturbation passes the 
on-ramp, it triggers ST which spreads upstream in the 
course of time, (v) Downstream, ST eventually relaxes 
to free traffic, (vi) ST often persists over several hours, 
(vii) The transition from ST to free traffic occurs at a 
lower density and higher average velocity than the in- 
verse transition (hysteresis effect). 

Property (i) is related to lane-changing and requires a 
multi-lane model for its description, e.g. Q. In order 



to explain the other characteristic properties of ST, we 
will propose a macroscopic, effective one-lane model that 
was derived from a gas-kinetic level of description and 
treats all lanes in an overall manner. The model is also 
in agreement with other empirical findings |^],[l2| like the 
existence of metastable states, the typical propagation 
velocity of upstream jam fronts (between —10 and —20 
kilometers per hour), and the characteristic outflow Q ut 
from traffic jams of 1600 up to 2100 vehicles per hour and 
lane (depending on the road and weather conditions, but 
not on the initial conditions or the surrounding traffic 
density) @. 

Our model is based on a kinetic equation for the phase- 
space density p(x,v,t), which corresponds to the spatial 
vehicle density p(x,t) times the distribution P(v;x,t) of 
vehicle velocities v at position x and time t ||. (For an 
introduction to gas- kinetic traffic models see Ref. jL4j.) 
The kinetic equation has some similarities to the gas- 
kinetic Boltzmann equation for one-dimensional dense 
gases with the vehicles playing the role of molecules. 
However, there are also some features specific to traf- 
fic. Drivers want to accelerate to their respective desired 
velocities giving rise to a relaxation term that violates 
conservation of momentum and kinetic energy. More- 
over, when approaching a slower car that cannot be im- 
mediately overtaken, one has to decelerate while the car 
in front remains unaffected. This leads to an anisotropic 
interaction. Finally, the reaction of the drivers depends 
on the traffic situation ahead of them, making the inter- 
action non-local. 

The model equations for the lane-averaged vehicle den- 
sity p(x,t) = Jdv p(x,v,t) and the average velocity 
V(x, t) = p^ 1 J dv vp(x, v, t) are 
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where we use the notation f a (x,t) = f(x a ,t) with / e 
{p, V, 9} and an advanced 'interaction point' x a specified 
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later. Without on- or off-ramps, the density equation (|l]) 
is just a one-dimensional continuity equation reflecting 
the conservation of the number of vehicles. Thus, the 
temporal change dp/dt of the vehicle density is just given 
by the negative gradient —dQ/dx of the lane-averaged 
traffic flow Q = pV. Along on-ramps (or off-ramps), the 
source term Q rmp / (nL) is given by the actually observed 
inflow Q rmp > from (or outflow Q rmp < to) the ramp, 
divided by the merging length L and by the number n 
of lanes. The inflow has an upper limit that depends on 
the downstream flow on the main road ||l5|] . 

The velocity equation (|^) contains the velocity vari- 
ance 9(x, t) = p^ 1 Jdv [v — V(x, t)] 2 p(x, v, t). Instead of 
deriving a dynamic equation for 9 from the kinetic equa- 
tions, we use the constitutive relation 9 = A(p)V 2 with 
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where A = 0.008, AA = 0.015, p c = 0.28p max , and 
Ap = 0.1p max These coefficients can be obtained 

from single-vehicle data. Unfortunately, no such data 
were available for the motorway considered in [0, but 
similar values were obtained for another motorway fl6}| . 

The first term on the rhs of Eq. (^) is the gradient of the 
'traffic pressure' p9. It describes the kinematic dispersion 
of the macroscopic velocity in inhomogeneous traffic as a 
consequence of the finite velocity variance. For example, 
the macroscopic velocity in front of a small vehicle clus- 
ter will increase even if no individual vehicle accelerates, 
because the faster cars will leave the cluster behind. The 
second term denotes the acceleration towards the (traffic- 
independent) average desired velocity Vq of the drivers 
with a relaxation time r G [10 s, 50 s]. Individual varia- 
tions of the desired velocity are accounted for by a finite 
velocity variance. The third term of the rhs of Eq. (j^) 
models braking in response to the traffic situation at the 
advanced 'interaction point' x a = x + 7(l/p m ax + TV). 
In dense traffic, where most drivers maintain the safety 
distance TV, this point is about 7 vehicle positions in 
front of the actual vehicle position a;. The average safe 
time headway T is of the order of one second. For the 
'anticipation factor' 7, we assume values between one 
and two. The braking deceleration increases coulomb- 
like with decreasing gap (l//O a — l//O m ax) to the car in 
front (l/p a being the average distance between succes- 
sive vehicle positions, l//9 max the average vehicle length, 
and p max the maximum density). In homogeneous dense 
traffic, the acceleration and braking terms compensate 
for each other at about the safe distance. In general, the 
deceleration tendency depends also on the velocity differ- 
ence to the traffic at the interaction point. A gas-kinetic 
derivation leads to the 'Boltzmann factor' lOl 
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where Sy = (V — V a )/\J9 + # a is the dimensionless ve- 
locity difference between the actual location x and the 
interaction point x a . In homogeneous traffic, we have 
-B(O) = 1. If the preceding cars are much slower (i.e. 
5y ^S> 0), it follows B(Sv) — 2o~y- In the opposite case 
(i.e. 6y "C 0), we have B(Sv) ~ 0. That is, since the 
distance is increasing, then, the vehicle will not brake, 
even if its headway is smaller than the safe distance. 

In contrast to previous approaches, the above macro- 
scopic traffic model explicitly contains an anisotropic, 
non-local interaction term B(Sy)- This is not only es- 
sential for a realistic treatment of situations with large 
gradients of p(x, t) or V{x, t), but also for an efficient and 
robust numerical integration. Moreover, the prefactor of 
B has now been obtained from the plausible assumption 
that, at high densities, the time headway between suc- 
cessive vehicles is T. Finally, all model parameters arc 
meaningful, measurable, and have the correct order of 
magnitude. 

Our simulations have been carried out with an explicit 
finite-difference integration scheme and the following pa- 
rameter values: Vq — 128km/h, p max = 160 vehicles/km, 
T = 1.6 s, t = 31s, and 7 = 1.0. The response of equi- 
librium traffic to localized disturbances is similar to the 
Kerner-Konhauser model |t]]. For densities p < p c \ and 
p > p c 4, homogeneous traffic is stable, and for a range 
Pc2 < P < Pc3 of intermediate densities, it is linearly un- 
stable, giving rise to cascades of traffic jams ('stop-and- 
go traffic'). For the two density regimes p c \ < p < p c i 
and p C 3 < p < p c 4 between the stable and the linearly 
unstable regions, it is metastable, i.e., it behaves nonlin- 
early unstable with respect to perturbations exceeding a 
certain critical amplitude, but otherwise stable. For the 
self-organized density pj am inside traffic jams we find a 
typical value p jam > p c4 @. 

Now, we will discuss synchronized flow. Figure 1 shows 
the simulation of freeway traffic near an on-ramp during a 
'rush-hour', where we assumed that the flow downstream 
of the on-ramp almost reaches the maximum equilibrium 
flow ('capacity limit') Q ma x- The upstream boundary 
condition at position xq = —6 km was specified in accor- 
dance with the equilibrium flow-density relation for free 
traffic (dotted lines in Fig. 3, before the maximum of the 
curve) with flows according to Fig. 1(c). We started with 
a high main flow that is monotonically decreasing in the 
course of time. At x = km, an on-ramp with merging 
length L — 300 m injects an additional time-dependent 
inflow Q rmp into the freeway. This on-ramp flow was 
assumed to have a short and tiny peak at t = lOmin 
[Fig. 1(c)]. As a result, a wave of denser traffic propa- 
gated downstream, thereby gaining a larger amplitude, 
and eventually propagated upstream again with about 
— llkm/h. Once the perturbation reached the ramp, 
dense traffic (of about 48 vehicles/km) with relatively 
high flows (1600 vehicles/h) corresponding to V =33 
km/h built up in the upstream direction. Although the 
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flow from the main road was gradually decreased for t > 
30 min, it took more than 100 additional minutes, until 
the region of congested traffic vanished. 

All these features agree with the experimental obser- 
vations of ST described in Ref. Q. There, a peak on the 
on-ramp flow was observed at about 7:15 am. The tran- 
sition to ST was first detected at 7:16 am as a short dip 
of the velocity 700m downstream from the on-ramp (de- 
tector D3 in H). At about 7:22 am, the front reached a 
detector (D2) 200m upstream of the ramp (corresponding 
to a mean propagation speed of — llkm/h), and propa- 
gated slower to the next detector Dl (700m upstream). 
While the perturbation at detector D3 lasted only a few 
minutes, it was followed by nearly 2 hours of congested 
traffic (V w 30 km/h, Q w 1500 vehicles/h) at the de- 
tectors D2 and Dl. 

It turned out that, apart from fluctuations, the simu- 
lated velocities and flows obtained at the detector posi- 
tions x = -0.7km (Dl), x = -0.2km (D2), x = 0.7km 
(D3), and x — 1.5 km (D4) (cf. Fig. 2) are in almost 
quantitative agreement with all features of ST as dis- 
played in Figs. 2(c), 2(b), 2(a), and 2(d) of Ref. [|. In 
particular, the model reproduces the drop of the velocity 
to about 30 km/h for up to two hours, while the flow 
is reduced by only 20%. Moreover, after the transition 
to free flow, the velocity is higher and the flow is lower 
than immediately before the transition to synchronized 
flow, both in the measurements and the simulation. Fi- 
nally, in Fig. 3 we depict the relaxation to free traffic 
downstream of the ramp by flow-density diagrams [see 
also Fig. 1(b)]. The results agree well with the empirical 
traffic data presented in Fig. 3(c) of Ref. [||. 

Our results suggest the following interpretation of the 
phase transition to ST. Initially, the homogeneous flow 
Qmain upstream of an on-ramp is stable, while the higher 
downstream flow Qdown = Qmain + Qxmp/ n is metastable 
(n = number of lanes). A perturbation of the ramp flow 
Qrmp triggers a stop-and-go wave, which travels down- 
stream as long as it is small and upstream as it becomes 
larger, as is known from 'localized clusters' M. Now, 
assume the downstream front of the cluster would pass 
the on-ramp. Then, since Qmain [Fig. 1(c)] is lower than 
the characteristic outflow Q ou t from a jam (being of the 
order of 2000 vehicles/km), the cluster would eventually 
vanish. However, during its lifetime, the cluster would 
continue to emit the flow Q ou t, leading downstream of 
the ramp to a flow Q ut + Qrmp/"- > Qmax- As a conse- 
quence, as soon as the perturbation reaches the on-ramp, 
it induces congested traffic with a standing downstream 
front just at the end of the ramp. With an observed out- 
flow Qout ^ Qout from ST jl7| , the average flow upstream 
is given by 

Qsync — Qout Qrmp/" ■ (5) 

Now, consider the density p S ync defined by Q sy nc = 
Qe(psync) in the congested part of the equilibrium flow- 



density relation Q c (p) (dotted lines in Fig. 3, behind the 
maximum of the curves) . If homogeneous traffic is (me- 
ta-)stable at /9 S ync, the on-ramp induces ST, otherwise 
it induces dynamically changing states. The restriction 
Qrmp < Qout/2 |ll]] (corresponding to every second vehi- 
cle on the right freeway lane stemming from the on-ramp) 
implies Q syn c > (1 - j^)Qout and p sync < p jam , so that 
synchronized flow is significantly higher than the flow in- 
side traffic jams. 

We have proposed a macroscopic traffic model based 
on a gas-kinetic level of description that allows to de- 
scribe the empirically observed features of traffic flows. 
This Letter focussed on the simulation and interpreta- 
tion of ST, which is probably the most common form 
of congested traffic. We have triggered ST by a small 
peak in the inflow from an on-ramp, when the down- 
stream flow was close to freeway capacity. Synchro- 
nized traffic eventually resolved in downstream direction, 
but spread in upstream direction. It persisted for more 
than one hour, although the main flow was steadily re- 
duced. We also performed simulations without peaks, 
leaving everything else unchanged. In these cases, we ob- 
tained free traffic flow. This confirms that the proposed 
model can describe the hysteretic and bistable proper- 
ties of real traffic. Our interpretation of ST underlines 
the crucial role of the characteristic outflow Q ou t from 
congested traffic for traffic dynamics. The simple cri- 
terion Q ou t + Qrmp/"- > Qmax for the formation of ST 
can be useful for determining bottlenecks of the exist- 
ing road infrastructure as well as for planning efficient 
freeway networks. 
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FIG. 2. Temporal evolution ol (a) the average velocity and 
(b) the traffic flow per lane at four cross sections of the free- 
way near the on-ramp. In front of the on-ramp (x < 0), ST 
exists for a certain time interval. Downstream (x > 0), the 
traffic situation recovers towards a freely flowing state. The 
simulated overshooting at the beginning of the breakdown of 
average velocity is in agreement with empirical observations 
(cf. Fig. 1(b) in Ref. §). 
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FIG. 1. Spatio-temporal evolution of the lane-averaged 
density after a small peak of inflow from the on-ramp. The 
on-ramp merges with the main road at x = km with a merg- 
ing length of 300m. Traffic flows from left to right. In (a), 
the parabolically shaped region of high density corresponds to 
ST. Plot (b) shows the formation of this state in more detail. 
The time-dependent inflows Q ma in at the upstream boundary 
and Qrmp/n at the on-ramp are displayed in (c). 
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FIG. 3. Traffic dynamics in the flow-density plane (a) 0.2 
km upstream of the on-ramp and (b) , (c) at two downstream 
cross sections. The solid lines with the symbols (□) corre- 
spond to the simulation results of Fig. 1. All the trajectories 
start at p — 17 vehicles/km and Q = 1770 vehicles/h. The 
dashed line represents the equilibrium relation Q e (p) of the 
model. The vertical dotted lines indicate the stability limits 
Pel, p C 2, Pc3, and p C 4 (determined numerically). 



4 



